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Abstract — An almost ubiquitous assumption made in the 
stochastic-analytic approach to study of the quality of user-service 
in cellular networks is Poisson distribution of base stations, often 
completed by some specific assumption regarding the distribution 
of the fading (e.g. Rayleigh). The former (Poisson) assumption is 
usually (vaguely) justified in the context of cellular networks, by 
various irregularities in the real placement of base stations, which 
ideally should form a lattice (e.g. hexagonal) pattern. In the first 
part of this paper we provide a different and rigorous argument 
justifying the Poisson assumption under sufficiently strong log- 
normal shadowing observed in the network, in the evaluation of a 
natural class of the typical-user service-characteristics (including 
path-loss, interference, signal-to-interference ratio, spectral effi- 
ciency). Namely, we present a Poisson-convergence result for a 
broad range of stationary (including lattice) networks subject to 
log-normal shadowing of increasing variance. We show also for 
the Poisson model that the distribution of all these typical-user 
service characteristics does not depend on the particular form 
of the additional fading distribution. Our approach involves a 
mapping of 2D network model to ID image of it "perceived" by 
the typical user. For this image we prove our Poisson convergence 
result and the invariance of the Poisson limit with respect to the 
distribution of the additional shadowing or fading. Moreover, in 
the second part of the paper we present some new results for 
Poisson model allowing one to calculate the distribution function 
of the SINR in its whole domain. We use them to study and 
optimize the mean energy efficiency in cellular networks. 

Index Terms — Wireless cellular networks, Poisson, Hexagonal, 
convergence, shadowing, fading, spectral/energy efficiency, opti- 
mization 

I. Introduction 

Cellular networks are being extensively deployed and up- 
graded in order to cope with the steady rise of user-traffic. This 
has created the need for new and robust analytic techniques 
to study the quality of user-service. The ability to tractably 
model and calculate the quality of user-service related to the 
signal-to-interference-and-noise ratio (SINR) will serve as the 
motivating force behind the work presented here. 

In order to derive analytic techniques, various mathematical 
models have been proposed. A common and simplifying model 
assumption is that the base stations are located according to a 
Poisson process in the plane. In the first section of this paper 
we recall such a model with shadowing and/or fading, and 
present a simple yet very useful result, upon which the bulk of 
the remaining work here hinges. This result involves mapping 
the point process on R? (that models the locations of base 
stations) through the distance-loss function and the shadowing 
and/or fading variables (that model the propagation losses 
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between these stations and a typical user), to a point process of 
propagation losses on R + , which are experienced by this user. 
It allows one to study all typical-user characteristics, which 
can be expressed in terms of its propagation losses (or received 
powers, e.g. SINR, etc). We observe that this new process also 
forms a Poisson process, regardless of the shadowing and/or 
fading distribution, and is characterized only by the moment 
of this distribution of order 2/(3, where is the distance-loss 
exponent. 

In the context of an actual deployment of cellular networks, 
lattice (e.g. hexagonal) models for the base station placement 
are usually thought of as more pertinent. However, perfect 
lattice models do not seem to allow analytic techniques for the 
study of the SINR-based characteristics. Hence, the Poisson 
model is used and justified by positioning "irregularities" of 
the network. It is also considered as a "worst-case" scenario 
due to its complete randomness property. Although the va- 
lidity of these arguments alone may be questioned, we seek 
to support the Poisson assumption with a powerfully new 
convergence result when there is sufficiently strong log-normal 
shadowing in the network. Namely, provided a network is 
represented by a sufficiently large homogeneous point pattern 
(whose definition will later be made precise, and includes 
general lattices and various "perturbed lattice" models), we 
show that as the variance of log-normal shadowing increases, 
the resulting propagation losses between the stations and the 
typical user form a stochastic process that converges to the 
aforementioned non-homogeneous Poisson process on R, due 
to the Poisson model. In other words, the actual (large but not 
necessarily Poisson) network is perceived by a typical user as 
an equivalent (infinite) Poisson network, provided shadowing 
is strong enough, of logarithmic standard deviation greater 
than approximately lOdB, as shown by numerical evidence. 
This is a realistic assumption for outdoor and indoor wireless 
communications in many urban scenarios. 

This result rigorously justifies using a Poisson point pro- 
cess to model typical-user characteristics in a wide range of 
network models, which includes the hexagonal model and a 
large class of perturbed lattice models, thus adding theoretical 
weight to the work being done under the Poisson assumption. 

After stating this convergence result, we derive important 
complementary analytic tools for the study of the distribution 
of the SINR in the Poisson model, which remain valid for any 
distribution of the shadowing and/or fading random variables. 
We use them to investigate the spectral and energy efficiency 
in cellular networks. In particular, we evaluate the mean 
energy efficiency as a function of the base station transmit 
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power, which allows one to optimally tune this latter power. 
Poisson and hexagonal networks with and without shadowing 
are compared in this context. 

Related work: It is beyond the scope of this short intro- 
duction to review all works with Poisson model of wireless 
networks. Some results and further references may be found 
in [1, 2], More specifically, our results presented in this paper 
allow to extend to a general fading distribution the explicit 
expressions for the distribution function of the SINR derived 
in [3], where Rayleigh fading is assumed. Moreover, using 
the Laplace/Fourier analysis, we manage to characterize (no 
longer explicitly) the SINR over its entire range, whereas 
the aforementioned explicit expressions are valid only for 
SINR > 1. The result of Lemma 1 appeared in [4]. The infinite 
Poisson model was statistically fitted in [5] to some real data 
provided by an operator regarding propagation losses in order 
to estimate the parameters of the propagation loss using a sim- 
ple linear regression model. The spectral and energy efficiency 
in hexagonal networks without shadowing was studied in [6]. 
Finally, the convergence result presented here is in the spirit of 
classical limit theorems of point processes, which are detailed 
in [7, Chapter 11]. In particular, these theorems show that 
under specific conditions, the repeated superposition, thinning 
or translation of points of a point process will result in the 
point process converging in the limit to a Poisson process. 

The remaining part of this paper is organized as follows. In 
section II we present the basic result on the network mapping 
to the process of propagation losses. The main convergence 
result is presented in Section III and its proof is given 
Appendix A. We revisit the Poisson model in order to study 
the SINR distribution in Section IV. 

II. Infinite Poisson model 

For motivation purposes, we first recall the usual cellular 
network model based on the Poisson process. In particular, we 
model the geographic locations of the base stations with an 
homogeneous Poisson point process $ = {Aj} ieN of intensity 
A on R 2 , which we refer to as the infinite Poisson model. 

We assume that the (typical) user is located at the origin 
without loss of generality due to the stationarity of {J^i} igN . 
Let I (Xi) be the distance loss between a base station at Xi 
and the user, where I (•) is given by 

l{x) = {K\x\f (1) 

for two given positive constants K and (3. When we incorpo- 
rate shadowing (and/ 'or fading), the propagation loss is 

L Xi - ~c 

where {S x } xeR2 is a collection of independent and identically 
distributed (iid) positive random variables. We will sometimes 
write also = Si to simplify the notation. Let 5 be a 
random variable having the same distribution as any S x . 

Note that the power received at the origin from the station 
Xi, transmitting with power Px t , is equal to 

P Xi = Px>S Xz 
Lx ( l(Xi) " 



In this paper we do not assume any power control, i.e., 
Pjq = P for some given positive constant P > 0. In this 
case (as well in a more general case of iid emitted powers), 
including Px t in the associated shadowing random variable, 
we retrieve an equivalent model in which the shadowing 
is Sxi '■= PxiSxi an d the transmitted powers Pxi = 1. 
Henceforth we assume, without loss of generality, that the 
transmitted powers are all equal to one, while keeping in mind 
that the shadowing random variables now include the effective 
transmitted powers. This transformation will slightly simplify 
our notation. However, when studying the energy efficiency 
in Section IV-B4 we will reintroduce emitted powers to our 
model. 

A. Mapping of the propagation losses of the typical user from 

m 2 to r+ 

Denote by Af = {Lxi} ie ^ the process of propagation losses 
experienced by the typical user with respect to the stations 
in $. We consider Af a point process on R + . Note that all 
characteristics of the typical user, which can be expressed 
in terms of its propagation losses (or received powers, under 
the aforementioned assumption on emitted powers, e.g. SIR, 
SINR, spectral and energy efficiency, etc) are determined by 
the distribution of Af. This motivates the following simple 
result that appeared, to the best of our knowledge for the first 
time, in [4]. In order to make this presentation more self- 
contained we present it with a proof. 

Lemma 1: Assume infinite Poisson model with distance- 
loss (1) and generic shadowing (and/or fading) variable 
satisfying 

E[Sf]<oo. (2) 

Then the process of propagation losses Af experienced by the 
typical user is a non-homogeneous Poisson point process on 
M. + with intensity measure 



where 



A ([0,i)) :=E[Af([0,t))]=ot* 



\icE[Sv] 



K 2 



(3) 



(4) 



Proof: The point process Af may be viewed as a trans- 
formation of the point process $ by the probability kernel 

p(x, A) —P € Aj , xeR 2 7 AeB(R+) 

By the displacement theorem [1, Theorem 1.10], the point 
process Af is Poisson on R + with intensity measure 



A([0,i)) = A 
= A 



IJx) 
S 



e [o, t) dx 



l(x) 



< t \ dxP s (ds) 



( AttE S? 



which completes the proof. 
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Remark 2: Note that the distribution of N is invariant with 
respect to the distribution of the shadowing/fading S having 
same given value of the moment ~E[S 2 ^]. This means that 
the infinite Poisson network with an arbitrary shadowing S is 
perceived at a given location statistically in the same manner 
as an "equivalent" infinite Poisson with "constant shadowing" 
equal to s const — (E[S* 2 ^])^/ 2 (to have the same moment 
of order 2//3). The model with such a "constant shadowing" 
boils down to the model without shadowing (S = 1) and the 
constant K replaced by K = K/y / 'E[S 2 ^]. 

The above result requires condition (2), which is satisfied 
by the usual models, as e.g. the log-normal shadowing or 
Rayleigh fading. Henceforth, when not specifying any par- 
ticular distribution of S we tacitly assume that condition (2) 
holds true. 

III. Convergence results under log-normal 

SHADOWING 

In this section we derive a powerful convergence result 
rigorously showing that the infinite Poisson model can be used 
to analyse the characteristics of the typical user in the context 
of any fixed ( deterministic!) placement of base stations, meet- 
ing some empirical homogeneity condition, provided there is 
sufficiently strong log-normal shadowing. 

A. Model description 

Let 4> = {Xi}^] be a locally finite deterministic point 



We consider also the analogous process of propagation losses 



pattern on 



For < A < oo, as r — > oo we require the 



empirical homogeneity condition 



A. 



(5) 



Note that the above condition is satisfied by any lattice (e.g. 
hexagonal) pattern <\>, as well as by almost any realization of 
an arbitrary ergodic point process. 

Let the shadowing between the station AQ and the 
origin be iid (across i) log-normal random variables 

Sf 5 - cxp(- ( r 2 /2 + aZi), (6) 

where Z{ are standard normal random variables. Note that for 
such S\ a) = S {a \ we have E[S] = 1 and E [(S^) 2 ^} = 
exp[a 2 (2-/3)//3 2 ]. 

Consider the distance-loss model (1) with the constant K 
replaced by the function of a 

a 2 {P-2 



K(&) = K exp 



2/3 2 



(7) 



where K > and /3 > 2. 

As in Section II, we consider the point process on R + of 
propagation losses experienced by the typical user with respect 
to the stations in <b 



M {0 



S 



X, G 



(8) 



This also means that the path-loss from a given station X expressed in dB, 
i.e., dB(Lx(y)), where dB (x) := 10 X log 10 (x) dB, is a Gaussian random 
variable with standard deviation ct c ib = oTO/loglO, called logarithmic 
standard deviation of the shadowing. 



#):=W,.< m < S .,,, £ 
I S- 



(9) 



where the stations in eb that are closer than a a and farther than 
b a are ignored, for all sequences < a a < b„ < oo satisfying 

log(a CT 



iog(M 



0, 



oo. 



(10) 



(11) 



B. Main result 

We present now our main convergence result. 

Theorem 3: Given homogeneity condition (5), then Af^ 
converges weakly as a — > oo to the Poisson point process 
on R + with the intensity measure A given by (3) with a = 
Xk/K 2 . Moreover, fif^ also converges weakly (a — > oo) 
to the Poisson point process with the same intensity measure, 
provided conditions (10) and (11) are satisfied. 

The proof of Theorem 3 is deferred to Appendix A. 

Remark 4: The above result, in conjunction with Lemma 1, 
says that infinite Poisson model can be used to approximate the 
characteristics of the typical user for a very general class of 
homogeneous pattens of base stations, including the standard 
hexagonal one. The second statement of this result says that 
this approximation remains valid for sufficiently large but finite 
patterns. 

Remark 5: The distance-loss model (1) suffers from having 
a singularity at the origin. This issue is often circumvented by 
some appropriate modification of the distance-loss function 
within a certain distance from the origin. The second statement 
of Theorem 3 with a a = const > shows that such a 
modification is not significant in the Poisson approximation. 

To illustrate Theorem 3 and obtain some insight into the 
speed of convergence we used Kolmogorov-Smirnov (K-S) 
test, cf [8], to compare the cumulative distribution functions 
(CDF) of the SIR of the typical user (see Section IV-A2) in 
the infinite Poisson model versus hexagonal one consisting of 
30 x 30 = 900 stations on a torus. We found that for 9/10 
realizations of the network shadowing the K-S test does not 
allow to distinguish the empirical (obtained from simulations) 
CDF of the SIR from the CDF of SIR evaluated in the infinite 
Poisson model with the critical p-value fixed to a = 10% 
provided adB ^ lOdB. On Figure 1 we present a few examples 
of these CDF's. Similar numerical study was done for the CDF 
of the path-loss with respect to the serving station in [5]. 

IV. Poisson model revisited 

We now return to the model outlined in Section II and 
investigate it further. We show that Lemma 1 is very useful 
in studying random quantities of the network, independently 
of the shadowing/fading distribution. We do not pretend to 
develop a complete theory, which is well beyond the scope 
of this paper. As a general remark let us emphasize only that 
many, already known, results (cf Related work in the Introduc- 
tion), were originally derived under specific assumptions on 
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Fig. 1. Empirical CDF of SIR simulated in hexagonal network with 
shadowing and their Poisson approximations. 



the distribution of shadowing and/or fading. By Lemma 1 they 
necessarily remain valid for a general distribution, with the 
appropriate specification of the value of the moment E[S 2 ^]. 

Our specific goal in this section is to present a new 
representation of distribution function of the SINR, which in 
consequence will allow us to study the spectral and energy 
efficiency — two important engineering characteristics of the 
network. To the best of our knowledge, they have not yet 
been studied in the Poisson model. Moreover, as another 
illustration of our convergence result, we will compare these 
characteristics in Poisson and hexagonal network, with this 
latter studied by simulations. 

Unless otherwise specified all results of this section regard 
Poisson model of Section II. 

A. Path-loss and SIR 

The novelty of our approach consists in representing the 
SINR in terms of the path-loss to the serving station and the 
respective SIR, rather than, as usual, the interference. 

1 ) Path loss: The weakest propagation loss denoted by 



L = inf Lx 



(12) 



Pl (ds) 



2a 



at? 



(It. 



(13) 



This corresponds to a Frechet distribution with shape param- 
eter i and scale parameter a? . 

2) Interference factor and SIR: The interference factor is 
defined by 

/ • 1 



L 



1 



where L is the path-loss factor. It may be interpreted as the 
interference to signal ratio; that is the inverse of the SIR. 
Introducing / = X)ieN j^r ( can t> e interpreted as the total 
received power) we can write / = LI — 1. Also, / = LI' 
where 

may be viewed as the interference. 
Define 



<P0 ( z ) 



+ z'7 1 



(14) 



where 7(0;, z) — J* t a ~ 1 e~ t dt is the lower incomplete 
gamma function. Another representation of the function fp{z), 
used in evaluation of the Laplace transform of the / (cf proof 
of Corollary 8) is given in Appendix B. 

Here is a key technical result for our approach. 

Proposition 7: The Laplace transform of the interference 
factor conditional to the path-loss factor is equal to 



E 



L 



a -a[v^(z)-l]s? 



Proof: We have 



E [e 



~ zf \L = si = Efe 



— zsl 



L = s]. 



(15) 



(16) 



is often called path-loss factor. It may be interpreted as 
the propagation loss with the serving base station (that is, 
the one with strongest received power). Note by (3) that 
the number of points of the process Af = {Lxi} ie fq is Thus 
almost surely finite in any finite interval. Hence, the above 
infimum is almost surely achieved for some base station; that is 
L = miriigN Lx t - Moreover, Lemma 1 allows us to conclude 
the following characterisation of the distribution of L, which 
follows immediately from the well known expression for void 
probabilities of Poisson process P{L > t} = P{Af([0,t)) = 
0} = cxp[-A([0,i))]. 

Corollary 6: The CDF of L is equal to V{L < t} 
1 — cxp[— at 2 IP\ where a is given by (4). Consequently, the 
probability density of L is given by 



Observe that /' is a shot-noise associated to the point process 
obtained from the Poisson point process Af of the propagation 
losses by suppressing its smallest (nearest to the origin) point, 
which is located precisely at L. By the well-known property 
of the Poisson process, given L = s, this new process is also 
Poisson on (s, 00) with intensity measure 



M ((s, f)) = A ([0, t)) - A ([0, s}), te (s, +00) 



2a 2 

M (dt) = A (dt) = —tfi~ 1 dt : on (s, +00) 

P 

and by [1, Proposition 1.5] 



(17) 



-zsl 



L 



exp 
exp 



2a 

'J 



(e * - l)M(dt) 



dt 



(18) 



Now we calculate the integral on the right-hand side of the 
above equation by making the change of variable u := ^f, 
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that is 



(1 - e-^)tTi~ l dt 



-du 



= (zs)% / (1 - e-^vri^du 
Jo 

= sf|[-l + ^(z)] 

where the third equality is obtained by integration by parts. 
Combining the above equation with Equations (16), (18) gives 
the final result. ■ 
The following results can be derived from Proposition 7, 
although we will not use them in the remaining part of the 
paper. 

Corollary 8: The joint distribution of L and f is charac- 
terized by 



E[l{L> u} e" zf ] 



1 



-au & (ff3 (z) 



z € 



The unconditional Laplace transform of the interference factor 
is 

1 



E [e-'f] 



z e 



<pp 0) ' 

Proof: For the fist statement we have 

/>oo 

E [1 {L > u} e- zf ] = E [e- zf \L = s] V L {ds) 

J U 



(19) 



1 2 On 3 



<P0 (z) 



where for the second equality we use (13) and (15). This 
completes the proof of the fist statement. For the second 
statement, conditioning on the value of the path-loss factor 
yields 

E[e- Z/ ] - / E[e-* f \L = s]P L (ds) 

JS.+ 

where (•) is the distribution of L given by (13). Using the 
above equation and (15), we ascertain 



li P 

POO 2 



E[e- Zf ]= I e -*[V^z)-L\sP — s f-± e -as» fa 



2a 



1 pOO rj 2 

fp(z) Jo P 
1 



where the third equality stems from (pp(z) ^ which follows 
from (40) in Appendix B and the assumption (3 > 2. ■ 



Remark 9: For t > 1 the (complementary) CDF of the SIR 
admits the following explicit expression 



P{SIR > t} = P{f < 1/t} = 



cW) 



(20) 



where 



C'(f3) 



2n 



£sin(27r//3) 



2r(2//3)r(l - 2/P) 



and T(z) = J °° e~ z z t ~ 1 dt it the complete Gamma function. 
It was proved in [3] assuming exponential distribution of S. 
But in the infinite Poisson model, SIR is invariant with respect 
to A, and K and consequently by Lemma 1 it is also invariant 
with respect to the value of EfS" 2 /^]. Hence the result remains 
valid for arbitrary distribution of S, provided E[5 2 /' 9 ] < oo. 
Other results of [3], including these involving the superposition 
of independent Poisson models (called if-tier cellular network 
model) can also be appropriately generalized using Lemma 1 . 

B. Distribution of SINR, spectral and energetic efficiency 

1) SINR: In this section we are primarily interested in the 
signal to interference and noise ratio 



SINR 



N - 



E 



_L 

1 

i6N L x - 



NL + f 



(21) 



where N is the noise power. We will show how CDF of the 
SINR can be evaluated, which opens a way for the study of 
functionals of SINR. 

2) Spectral efficiency: An important characteristic of a 
wireless cellular network is its spectral efficiency, defined in 
the simplest case of additive white Gaussian noise (AWGN) 
and the optimal theoretical link performance as 

S := log (1 + SINR) . 

It tells us how many bits per second and per Hertz can be sent 
to the typical user of the network. 

3) Energy efficiency: Up to now we have considered a unit 
transmitted power. Assume now that base stations transmit 
some power P > 0. In fact, in order for a base station to 
be able to transmit this power to the mobile, it needs to be 
powered (i.e., consumes energy per second) at the level P' > 
P. Following [9], assume that these two quantities are related 
through a simple linear relation P' = cP+d for some positive 
constants c and d. The energy efficiency is defined by 



£ :=£{P) 



W log 1 



NL/P+f 



cP + d 



where W is the bandwidth expressed in Hz. Thus £ is equal to 
Wnog(l + SINR(P))/P', where SINR(P) takes into account 
the transmitted power, and P' is the consumed power. It tells 
us how many bits per second per Watt of consumed power can 
be sent by a base station to the typical user. Note that £(0) = 
f (oo) = and thus £(P) admits a non-trivial optimization 
in P. 
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4) Evaluation of the CDF of the SINR: Proposition 7 in 
conjunction with Corollary 6 completely characterizes the joint 
distribution of L and /. It does not, however, allow for an 
explicit expression for the CDF of the SINR in the whole 
domain (cf Remark 9). In this section we describe a practical 
way for numerical computation of this CDF. We will use it 
to study the s pectral and energy efficiency. In fact, we will 
compute the CDF of the random variable Y := NL + f, which 
is sufficient, in view of (21). 

Proposition 10: The cumulative distribution function of Y 
is given by 



p oo 

P(Y<x)= F s (x - Ns)P L (ds) 
Jo 



(22) 



where P l (da) 
P(f<y\L = s) 



is given by (13), 
may be expressed by 



and F s (y) 



F s (y) = l — / K(£ F h + iu)) cos ut du, (23) 

n Jo 

where 7 > is an arbitrary constant, 1Z is the real part if it 
complex argument and 



£f s (z) = 



= I _ I e -«# fl (*)-i]J 



(24) 



with (pp (•) is given by (14). 

Proof: The expression (22) is trivial. The expression (23) 
of the (conditional) CDF F s is based on the Bromwich contour 
inversion integral of the Laplace transform Cp s of 1 — F s . The 
expression (24) for this latter follows from (15). ■ 

Remark 11: As shown in [10], the integral in (23) can 
be numerically evaluated using the trapezoidal rule, with the 
parameter 7 allowing control of the approximation error. The 
function F s may be also retrieved from (15) using other 
inversion techniques. 

Remark 12: For t > 1 the (complementary) CDF of the 
SINR admits the following expression 



P{SINR > t} 



r(i + f) 



-r ! r(l-2/ffl-Ma-W 2 r' i 



dr 



(25) 

with a given by (4). It follows from [3, Theorem 1] (where the 
exponential distribution is assumed for S) and our Lemma 1. 
Note on Figure 2 that the expression is not valid for t < 1. 

5) Numerical examples: In what follows we assume the 
following parameter values. If not otherwise specified we take 
log-normal shadowing with logarithmic standard deviation 
<7dB = 12dB (cf the footnote on page 3). The parameters 
of the distance-loss model are K — 4250Km~ 1 , /3 = 3.52 
(which corresponds to the COST-Hata model for urban en- 
vironment with the BS height 30m, mobile height 1.5m, 
carrier frequency 805MHz and penetration loss 19dB). For the 
hexagonal network model simulations, we consider 900 base 
stations (30x30) on a torus, with the cell radius R = 0.26Km 
(i.e., the surface of the hexagonal cell is equal to this of the 
disk of radius R). System bandwidth is W = 10MHz, noise 
power TV = — 93dBm. 
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Fig. 2. CDF of the SINR for (finite) hexagonal network with and without 
(°dB = OdB) shadowing, finite Poisson as well as for the infinite Poisson 
model. The last in the legend curve corresponds to (20). 



120000 
100000 














hexsgona |F :r]l 
exagonal. n dB -i2dE 
Infinite Po ssor 
























80000 
60000 
40000 






































20000 









































Transmitted power [dBrn 



Fig. 3. Mean energy efficiency of finite hexagonal model with and without 
(°dB = OdB) shadowing as well as the infinite Poisson in function of the 
transmitted power. 



Figure 2 shows the CDF of the SINR assuming the base 
station power P = 58.5dBm. We obtain results by simulations 
for finite hexagonal network on the torus with and without 
shadowing as well as the finite Poisson network on the same 
torus. For the infinite Poisson model we use our method based 
on the inversion of the Laplace transform. For comparison, we 
plot also the curve corresponding the explicit expression (20), 
which is valid only for SINR > 1. 

Figure 3 shows the expected energy efficiency E(£) of the 
finite hexagonal model with and without shadowing as well 
as the infinite Poisson, assuming the affine relation between 
consumed and emitted power with constants c = 21,45 and 
d = 354.44W. 

On both figures we see that the infinite Poisson model gives 
a reasonable approximation of the (finite) hexagonal network 
provided the shadowing is high enough. In particular, the value 
of the transmitted power at which the hexagonal network with 
the shadowing attains the maximal expected energy efficiency 
is very well predicted by the Poisson model. 

Conclusions 

We present a powerful mathematical convergence result 
rigorously justifying the fact that any actual (including regular 
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hexagonal) network is perceived by a typical user as an 
equivalent (infinite) Poisson network, provided log-normal 
shadowing variance is sufficiently high. Good approximations 
are obtained for logarithmic standard deviation of the shad- 
owing greater than approximately lOdB, which is a realistic 
assumption in many urban scenarios. Moreover, the equivalent 
infinite Poisson representation is invariant with respect to 
an additional fading distribution. Using this representation 
we study the distribution of the SINR of the typical user. 
In particular, we evaluate and optimize the mean energy 
efficiency as a function of the base station transmit power. 

Appendix 

A. Proof of Theorem 3 

In order to simplify the notation we set n := a 2 . Moreover, 
without loss of generality we assume that n takes positive 
integer values. Also, it is more convenient to study the point 
process of propagation losses on the logarithmic scale, which 
we do in what follows. In this regard denote by Ai og the image 
of the measure A given by (3) with a = Xtt/K 2 through the 
logarithmic mapping 

A log ((-oc,s]) := l(log(t) < s)k(dt) = -^ef 

for sel. 

For a given i, we first observe by (6) that 



u„(s,\Xi\) :=P 

= P 

= G 



log 



Z, < 



K{nf\X t \P 



S 



(«) 



< s 



s-piog^X^-n/p 



S -/31og(2rpQ|)-n//3 



(26) 



where G is the CDF of the standard Gaussian random variable. 

Let B (r) = {x E W 2 : \x\ < r}. We now need to derive 
two results. 

Lemma 13: 

lim / v n {si \x\)dx = lim / v n (s,\x\)dx 

n ^°° J B (b„)\B (a„) n ^ co iE2 

= A log ((-oo,s]) 

provided that a n and b n satisfy (10) and (11). 

Proof: We start with the integral over IR 2 in polar form 



/ 



v n (s, \x\)dx = 2tt / rG 
o 



s-/3log(Kr) -n//3 



dr 



Introduce a change of variables 

s — f3 log(ATr) — nj (3 



t 



1 

K 



exp 



(s — t\fn — n/P^j 



dt 
dr - 



P dr 
\fn r 
Jnr 



P 



dt. 



hence 

i' oo 

2tt / rG 
Jo 



s — P \og(Kr) — n/P 



dr 



= 2tt 



= 2tt 



OC 



K 2 



exp 



P 



(s — t\fn — nj /?) 



G(t)^dt 



j(s-n/P) 
K 2 



exp 



P 



G(t)dt 



Noting that G(~t) = 1 - G(t) := G(t), apply integration by 
parts on the above integral 



exp 



-2t^ph 
~P~~ 



G(t)dt = I exp 

— oo 

p 



2y/7l 



exp 



2ty/n~ 
2ty/n 



+■ 



P 



pOO 




/ exp 

' —oo 


P 



G(t)dt 

OO 

G(t) 

— OO 

g(t)dt, 



where g(t) is the normal probability density, and, via the 
inequality 

-y 2 

e- x2 dx<-^-—, y>0, (27) 

>y (i + y) 

(cf [11, Section 7.8]) the final integrated term on the left 
vanishes. Hence 



exp 



-2ty/n' 



2VTi J_ 



2y/n 



— exp 



2^/n 
which gives 



exp 



P 

exp 

j 

2n 
J 2 

2n 

J 2 



G(t)dt 

t 2 2ty/n 
~2 + ~_ 



dt 



2?r 



exp 



1 / 2JTi 



dt 



2n 



(28) 



2tt / rG 
>o 



s — p\og(Kr) — n/P 



7T 2s 



We proceed similarly with the other integral in Lemma 13 

v n (s, \x\)dx 



B (bn)\B {a n ) 

2tt I rG 



s- p\og(Kr)-n/P 



dr. 



The change of variables t = s 13 l ° e ^T^ gives 



r b„ 

2tt I rG 



s — P \og(Kr) — n/P 



dr 



= 2tt / exp 



P 



(s — t\fn — nj P) 



G(t)^dt 



2n- 



■ exp 



P 



$(s-n/P) 
K 2 



exp 



-2ty/n 
~1~ 



G(t)dt 
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where 



Moreover 



exp 



exp 



J_ 

J_ 
2^/n 



exp 



exp 



s-piog(Ka n )-n/p 
s/n 

i- p\og{Kb n )-n/f3 
\fn 

G(t)dt 
G(t)dt 
G{t) 



for some k > 0, whose value will be fixed later on. In the 
limit of n —> oo, the first summation in (31) disappears 



£ Ms,\x t \) 

Xie4>r\B (r kg ) 



2V" 





'2ty/n~ 


/ exp 

/ — U n 


P \ 



g(t)dt 



£ p 

= £ g 

Xi£<f>nB (r ko ) 

< <p(B (r ko ))G 
-> (n ->• c 



Z < 



s-^log^l^D-n/^ 



s-/31og(X|X 4 |)-n//3 



/n g(iT|X*|)-n//3 



s-/31og(/f|X|)-n//3 



G(t) 



2^h 



exp 



"2?i" 






/ exp 

J —u n 


; the 


conditions 



2^ 



where X* gives the maximum of G 
X £ 4>r\Bo(rk ) which exists since <J> is (by our assumption) a 
locally finite point measure. For the second summation in (31) 
we write u n (s, \Xi\) = u n (s, \x\^ffi), hence 



v n {s, \Xi\) 
1 



that the above integral converges and the integrated term 
disappears. The latter can be achieved, given inequality (27), 
in the limit as — u n and — v n both approach infinity, or equiv- Then the bounds 



1-4 



fcl J A 



v n {s, \x\^j-)dx. 
\x\ 



alently ^M|M + ^^00 and ^ s {Ka n) + ^ ^ 
as n — > 00, which agrees with conditions (10) and (11). 
Further conditions are revealed after the change of variable 
w = t — 2y/n/f3, yielding the integral 



Tk < \Xi\ < r k+1 



Tk+l 



and form of v n , which implies u n (s, \x\e e ) — v n (s — f3e, \x\ 
lead to the lower bound 



dt 



2tt 



f> _v2_ dw 
e 2 



2tt 



E E "«(«,l*il)> £ 
k=k Xi£ij>nA k 



k = k() 



whose limits of integration imply, in light of the earlier integral 

(28), that as n -> 00 the following /31 ° s ^ b " ) > 1//3 and and the upper bound 
l3\og(Ka n ) <- jy^g are required, of which both conditions (10) 00 
and (11) satisfy. 



E £ !/„(*, M < £ 

k=k Xi£4>nA k 



k—kn 



4>{A k 
\M 



<HAk 

\M 



Lemma 14: Assume (5), (10) and (11), then 



v n {s—fie, \x\)dx, 
(32) 

(33) 



lim 

n— yoo 



£ i^ii) 

X 1 e0n(S o (b„)\So(a„)) 

= lim } i/ n (s, \Xi\) = Aiog((-oo, s]). 



(29) 
(30) 



Moreover, we write 

0(A fc ) 0(S o (r fe+1 )) 



b(r*)) 



|4k 



|B (r fc+ i)|-|B (r fc )| 

0(Bo(r fc + i)) _ 0(B o (r fc )) \B (r k )\ 
|B (r fc+ i)| |Bo(r fc )| |B (r fe+ i)| 



1 



Proof: For fe > and a fixed e > 0, let r k = e ek and 
j4/c = Bo(rk+i) \ Bo(rk), and write the summation in (30) as 



\Bo(r k )\ 
\B (r k+1 )\ 
c/>(B (r k + 1 )) _ 0(B o (r fc ))^- 2e 



|Bo(r fc+ i)| 



|Bo(r*)| 



1-e 



-2e 



E ^(s,|Xi|)+£ E "n(s,\X t 



(31) 



and requirement (5) yields limfc^oo ^p^p = A. Hence, for 
any fixed S > 0, there exists a ko(S) such that for all k > /co, 
the bounds 

(l- j )A<^<(l + ^, 
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hold. Lower bound (32) becomes 

oo OO „ 

E E u ^ s ' I^D > E ( X ~ ^ A / ^(s - \x\)dx, 



k=k Q Xi£<pnA k 



k—kn 



= (1-S)X v n (s - 0e,\x\)dx 

J\x\>r k() 

Finally, Lemma 13 allows us to set a n = r/, g and b n = oo, 
hence 

oo » 

J^E E "MXiDZd-Qj^e 2 ^, 

k=k Xi£ipnA k 

and similarly the upper bound (33) becomes 

oo , 

n ^E E M S ,\X,\)<(l-S)^- 2 e 2 -^, 
k=k Xie4>nA k 

and indeed e —> and <5 ^ completes the proof of (30). 
The other result, (29), can be proved by a straightforward 
modification of the above arguments. ■ 
Proof of Theorem 3: We use a classical convergence 
result [7, Theorem 11.22V] which in our setting requires 
verification of the following two conditions (cf [7, (11.4.2) 
and (11.4.3)]) 

sup//;(A) ^ (n->oo). (34) 



and 



J2K(^)^Mo S (A) (n->oo), (35) 



for all bounded Borel sets A C K, where v l n {-) is the 
(probability) measure on R defined by setting v l n [(— oo, s]) := 
v n (s, \Xi\). The first condition, (34), clearly holds by (26) for 
any locally finite <fi. The second condition, (35) follows from 
Lemma 13 and 14, which establish the required convergence 
for A — (—oo, s] and any s £ K. This is enough to conclude 
the convergence for all bounded Borel sets. ■ 

B. Representation of the function tpp(z) given by (14) 

Lemma 15: The function tpp(z) given by (14) can be 
written as 



where 



r,1 -iM-!" 



7* (a, z) = ——j(a,z) 
T(a) 



E 

k=0 



r{a + k + l) 



is an analytic function on C x C. Moreover 

<pp(z)>0, V^>2,Vzet 



(36) 
(37) 

(38) 
(39) 

(40) 



function 7* defined by (38) (which is holomorphic on C x C; 
see [12, 6.5.4]) we obtain (37). The expansion (39) is given 
in [12, 6.5.29]. This expansion shows that if a > — 1 then for 
any z G R + , 7* (a, z) > 0, and therefore by (37) <pp (z) > 
for any /3 > 2,z e E+. ■ 
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Proof: Using (14) and recurrence formula [12, 6.5.22] 
we obtain (36). Introducing the modified incomplete gamma 



